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INFLUENCE OF FRICTION FORCES ON THE MOTION OF VTOL AIRCRAFT 


DURING LANDING OPERATIONS ON SHIPS AT SEA 
James C. Howard and David 0. Chin 
Ames Research Center 


SUMMARY 


The equations describing the friction forces generated during landing 
operations on ships at sea have been formulated. These forces depend on the 
platform reaction and the coefficient of friction. The platform reaction 
depends on the relative sink rate and the shock-absorbing capability of the 
landing gear. The coefficient of friction is assumed to vary with the surface 
condition of the landing platform and the angle of yaw of the aircraft relative 
to the landing platform. The landings contemplated are by aircraft of the 
VTOL type equipped with conventional oleo-pneumatic landing gears. Because 
aircraft of this type land with little or no forward speed, splnup forces have 
been neglected. Simplifications have been introduced to reduce the complexity 
of the mathematical description of the tire and shock-strut characteristics. 

It has been shown elsewhere that, for normal impact without tire bottoming, 
reasonable variations in the force-deflection characteristics of the tire have 
only a relatively small effect on the calculated behavior of the landing gear. 
Approximating the actual complicated force-deflection characteristics of the 
tire by a linear relationship appears to be adequate for practical purposes. 
Although the Internal friction forces in the shock strut have been Included in 
the landing gear model, experimental data suggest that these forces can often 
be neglected without impairing the validity of the model equations. By includ- 
ing only those characteristics of the tire and shock strut that contribute 
significantly to the generation of landing gear forces, a set of relatively 
simple equations is obtained. Nevertheless, the equations are considered 
adequate for practical purposes. 


INTRODUCTION 


Landing VTOL aircraft on ships at sea is complicated by the motion and 
surface condition of the landing platform, which moves in response to the 
motions of the sea, and which is exposed to a variety of atmospheric condi- 
tions. An understanding of the landing phenomena requires that a mathematical 
model, which Includes platform reactions and friction forces, be formulated 
and used to simulate the lending maneuver. In assessing the platform reaction 
forces likely to be encountered, the actual complicated force-deflection 
characteristics of the tire have been approximated by a linear relationship. 
Moreover, in the event that tire bottoming occurs, a linear segment approxima- 
tion which takes into account the Increased stiffness of the tire that -results 
from bottoming yields good results (Mllwltsky and Cook, 1953). A limited 
amount of experimental data obtained during drop tests at the Langley Research 



Center indicates that the behavior of the landing gear is relatively insensi- 
tive to variations in the air compression process , and that variations between 
isothermal and near adiabatic compression have only a secondary effect on the 
calculated behavior of the landing gear. Consequently, variations in the 
polytropic exponent have been neglected. The importance of friction forces 
derives from the fact that, if friction is not sufficient to prevent sliding, 
the aircraft may be damaged by colliding with adjacent structures. It is 
unlikely that VTOL aircraft landings will be attempted when the decks are 
awash, but during rough weather the landing platform will be exposed to rain 
and possibly sea spray which could reduce the coefficients of friction below 
safe levels. Another possibility to be considered is the likelihood that VTOL 
aircraft may produce a shower of spray as the jets Interact with the surface 
of the sea during the approach. If the resulting spray were to be blown over 
the landing platform, the friction forces w^ild be reduced significantly. In 
view of these considerations, it is important to determine the reactions and 
friction forces likely to be encountered during a variety of sea states and 
landing conditions. By including these forces in the equations of the mathe- 
matical model, the motion of the aircraft during the landing maneuver can be 
computed and the possibility of sliding determined. In formulating the fric- 
tion force equations, bending will be neglected and it will be assumed that 
each landing gear has a fixed orientation relative to aircraft body axes. In 
many cases, the influence of elasticity may be neglected without serious error, 
but in some Instances, particularly when the landing gear attachment points 
experience large displacements relative to the nodal points of the flexible 
system, the interaction between the deformations of the structure and the 
landing gear may be required to represent the system adequately. Expressions 
for the hydraulic, air-compression, and internal friction forces generated in 
the shock strut are derived in appendix B. 


KINEMATICS AND DYNAMICS OF AIRCRAFT-SHIP INTERACTIONS 


Forces and Moments 

Friction forces and platform reactions- Relative to the ship's coordinate 
system (XjXjX,), which originates at the center of gravity, the landing plat- 
form lies in the plane (X s - C) ■ 0, where C is negative constant. The 
X x X] plane of the ship's coordinate system corresponds to the plane to sym- 
metry of the vessel and the X 2 axis is normal to it. A platform coordinate 
system (xfxfxf) having axes parallel to the axes of the ship's coordinate 
system is used to specify vector components relative to the platform (see 
fig. 1). 


Subsequent to the instant at which wheel w makes contact with the land- 
ing platform, the aircraft is subjected to a reaction force Fjf which is 
normal to the landing platform and a friction force which is coplanar with the 
platform. Relative to platform axes, the friction force has components F* 
and F^. Hence, subsequent to wheel contact, the force vector Y* generated 
by wheel w has components 


F W - F?1 + Fj 3 + ijt 


( 1 ) 


2 



where I, J, and ft are a triad of 
mutually orthogonal unit vectors in 
the directions of the platform axes 
and the superscript w denotes the 
wheel being considered. 

Therefore) if the dynamic 
response of the aircraft subsequent 
to the wheel contact is being com- 
puted, the wheel reactions and fric- 
tion forces of equation (1) must be 
added to the inertia, thrust, aero- 
dynamic, gravitational, and shock- 
strut forces to complete the mathe- 
matical description of the force 
system. 

The corresponding wheel forces 
relative to aircraft body axes are 
fjJ where 

f i " d ij F j * i *J " 1 ' 2 * 3 < 2 ) 

and djj are direction cosines 
(appendix A). The force generated 
by each wheel is obtained by assign- 
ing the appropriate wheel number. 

In this and subsequent equations, 
the summation convention is <^sumed; 
that is, if in any term an index 
occurs twice, the term is to be 
summed with respect to that index for 
all admissible values of the index. 

Moments produced by friction 
forces and platform reactions- The 
ith moment component produced by the 
platform reaction and the friction 
force at wheel w is the vector 
cross-product of the position vector 
of wheel w and the friction force 
vector: 



Figure 1.- Systems of reference axes, 
including ship and platform. 


< - <*: - «> 

where x* are the aircraft body axes positional components of wheel w. It 
should be noted that 1, j, and k must be in cyclic order in equation (3). 
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From equation v.2) 


VS 


f" - d, F w 
J jn n 


Substitution of these values in equation (3) yields the following moment 
equation: 

•? • - *vV F « < 4 > 

Therefore, the individual components are 

M* “ < x X„ ~ x X > fW 

1 2 an J 2n n 

U W . W, w. v„W 

“ ^ x s d xn ~ x x d sn^ F n 

M w - (x w d - x w d )F W 
3 vx x 2 n *2 in' n 

n ■ 1,2,3 

When these equations are summed on n, the moment components assume the form 

+ d 32 F r + + d 22 F I * d 2 S F pl (5) 

^ + d xx F i> - << d ,x F x + d , a p : + ^ <*> 

< “ 1*>2X F T * d 22 F I + d „*T> ' X I< d XX F ? + d X2^ + d X, F T» < 7 > 


Orientation of Friction Force Vector 

Relative velocities- To determine the orientation of the friction force 
vector In the plane of the landing platform. It la necessary to know the 
relative velocity of each landing wheel at the Instant of touchdown. The 
friction force at each wheel will be in a direction opposite to the direction 
of the aircraft velocity relative to the landing platform at the point of 
wheel contact (Ijff, 19^). The point at which wheel w makes contact with 
the landing platform has a position vector R v relative to the origin of the 
ship's coordinate system. The angle In the plane of the landing platform 
between the direction of the axis and the velocity vector at wheel w Is 
denoted by k w (see fig. 2). When the magnitude of the wheel velocity rela- 
tive to the landing platform Is known, the angle te w is obtained as follows. 
The velocity of wheel w relative to the landing platform Is V w , where 

V* - v?i + V?J + VjK (8) 
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Figure 2.- Systems of platform reference axes. 


This velocity is obtained by subtracting the velocity of the point ot wheel 
contact on the platfora froa the wheel's velocity relative to the platform: 


fr' 

T ^ 


v“ - D(5 ct ♦ =.»*) - <? lhip ♦ :, hlp «i») 

In aatrlx notation these equations assuae the fora 

di2 




U + QX| - mV 


-I v ♦ rxY - pxY 


(9) 


W + PX? - qxY 


where U,V,W and P.Q.R are the components of the ship's linear and angular 
velocity vectors, respectively. 

For experiments of this type, the most important components of ship 
notion are heaving, pitching, and rolling. Ranee, by assuming that 


0 - V - 1 • 0 

the amount of computation la reduced, and equation (9) aaa 
form 




the simpler 




-w? 


( 10 ) 





Moreover , 



((? + wxr w ) ■ 

c 8 


r u + qxY * rx!T 


u u 

v + rxi - px3 


v w + p^2 - 


(ii) 


where V c _ is the velocity of the aircraft's center of gravity and r w is 
the position vector of wheel w relative to the origin of the aircraft body 
axes system (Etkln, 1972). Hence, in terms of the velocity components and 
the angle k w is (see fig. 2) 


k w - tan" 1 (vY/V*) (12) 

Equation (10) can be solved for the velocity components V* when the 
X* coordinates of the point of wheel contact on the platform are known. To 
determine these coordinates, it is first necessary to compute the position 
vector R C R where 


R cg - Xj 8 i ♦ X 2 g J + X® 8 '* 


(13) 


and X* 8 constitutes the locus of points on the trajectory of the aircraft 
center of gravity projected onto the plane (Xj - C) ■ 0, which is the plane of 
the landing platform. Hence the distance of any point on the projected tra- 
jectory from the origin of the ship's coordinate system is 

R cg - X® 8 i ♦ xJ 8 J + CK (14) 


where C is a negative constant. By substituting X* “ xj 8 and xi ■ 0 in 
equations (10) and (11), respectively, the velocity components of the air- 
craft's center of gravity relative to the landing olatform are obtained: 



Therefore, the velocity components V x 8 and V 2 8 are 


(15) 


V ? 8 - (d xl u ♦ d Jx v ♦ d Jx w - QC) 


(16) 


Vf 8 - (d l2 u ♦ dj jV + d, jW + PC) 


(17) 
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The velocities V^8 and Vf8 can be Integrated to yield the coordinates X*® 
and x£ Therefore, the required coordinates are 

x 5 8 • x 5f + / v i 8 dt (18) 

X 2 8 * X 29 + / V 2 8 dt 09) 

where Xj£ and X 2 f are Initial coordinates. 

The velocity of the aircraft center of gravity normal to the landing 
platfora is given by 

v' 8 - Ud ljU + d 2J v ♦ d,,v - (W + PXf 8 - QX* 8 )J (20) 

Since X 2 8 and X 2 8 are known from equations (18) and (19), respectively, 
equation (20) can be integrated to yield 



where Xjf 1^ the initial distance from the X,X 2 plane of the ship's coor- 
dinate system. The Initial values Xff, Xff, and Xff thip axes 

components. 

When initial conditions are given In Earth-fixed axes, a transformation 
from these axes to ship axes Is required. Given that X^® are Initial valuer 
In Earth-fixed axes, the corresponding ship axes components are (from 
appendix A) 



( 22 ) 


Vhe aquations (18) through (21) are modified in accordance vith this trsna 
formation, the coordinates of the aircraft trajectory assume the form: 



Let the origin locatior. of the landing platform with respect to the ship's 
coordinate system be X 1 * ■ ( X x o ’ x f 0 , C); then the height of the aircraft cen 
ter of gravity above tie landing platform is xf, where 

X* - (Xj® - C) 


Moreover, 

X l “ (Xx 8 - X* 0 ) 

xf - (X ® 8 - xf 0 ) 

Wheel w will make contact with the landing platform "hen (Howard, 1977) 

” - (d X3 X I + d 23 X I + d 33 X I ) (2 -> 

(see fig. 3). 



Figure 3.- Aircraft attitude relative to landing platform. 




When this condition is satisfied, the computed values of X x ^ and 
may be used to determine the point of contact of wheel w from 

X T “ X 5 8 + d ll*T + d 21 X 2 + d 31 X I (25) 

- X 2 8 + d X2 xY + d 22 x 2 + d J2 x* (26) 

The relative velocity components can now be evaluated by substituting from 
equations (25) and (26) in equation (10), and the angle k w determined from 
equation (12). 

Kinetic friction force- The friction fotce vector f” at wheel v is 
opposite the resultant velocity vector and is equal to the product of the 
coefficient of friction and the platform reaction. It is assumed that the 
coefficient of friction is a function of the angle t w , which Js the angle 
between the resultant velocity vector at wheel w and the line of intersection 
of the rolling plane of the wheel with the landing platform. Denoting the 
direction of this line by the unit vector p and the perpendicular direction 
by the unit vector n, the velocity components relative to these directions 
are (see sketch (a)) 


\ 



Vj| - vj •in<K* - 6^) - Vj sin t W (28) 
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where tc w is defined by equation (12) and 0^ is the aircraft angle of yaw 
with respect to platform axes. The velocity components at wheel w relative 
to platform axes are and and the resultant velocity in the plane of 
the platform is Vg. 

The equation for angle t w may be written in the following alternative 
forms (see eq. (A19)): 



Relative to these directions, the kinetics friction force vector assumes 

the form 


=w w , w~ . . w«,_w 

* p (cos t p + sin v u)F 3 


(29) 


Relative to aircraft body axes, the kinetic friction force and the platform 
reaction have components f^g, where 



(30) 


and 



cos 0, 


sin 0 3 sin e 2 cos 0 2 


^cos sin 0 2 -sin 0 2 


-sin 0 2 
sin 0 X cos 0 
cos 0 2 cos 0 



Kinetic friction moments- When the components of the force vector from 
equation (30) are substituted in equations (5) through (7), and the coeffi- 
cients Dij substituted for djj, the moment components are obtained as 
functions of the platform reactions. These are 
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•„V» '.itf** 


! 


I 



i 


M T “ I X 2< D 31^ W COS T * + D 32 uW 8in tW + D 33> 

- x«(D 21 u w cos t W + D 22 y W sin T W + D 23 )]F 3 (31) 

m" - [x"(D ii p W cos t w + D i2 u w 8in t W + D ia > 

- x"(D 31 y w cos t W + D J2 y W sin t W + D 33 )]F^ (32) 

M T " I*! (D 21 liW C ° 8 tW + D 22 WW 8i " T " + D 33 > 

- x 2 (D 31 p w cos t w + D 12 y w sin t W + D 13 )]F 3 (33) 

Experimental evidence indicates that, in the case of normal impact without 
tire bottoming, the assumption of a linear force-deflection relationship for 
the tire is adequate for practical purposes (Milwitzky and Cook, 1953). Hence, 
by assuming that each tire is a linear spring with an effective stiffness K e , 
the force resulting from a deflection 6 W can be expressed in the simple form 

F* - -K*6 W (34) 


Moreover, in the event that tire bottoming occurs, a linear segment 
approximation which takes into account the increased stiffness of the tire 
that results from bottoming yields good results. 


Subsequent to tire and shock-strut deflection, the moment arms become 
functions of the tire and shock-strut deflections; that is 


W , W «W Wv , /.W_W _ ,Wv \ 

x i " < x io “ * * ) + (* R " D 13 fi ) 

* 2 - (x* 0 - m w s w ) + (m w R w - D 23 6 W ) 

X? - <X? 0 - n w s w ) + (n w R w - D J3 « W ) , 


(35) 


where x” 0 are the coordinates of the axle at wheel w when the shock strut 
is in the fully extended position; l V , m w , and n w are direction cosines; s w 
is the shock-strut deflection, and R w is the radius of wheel w (see fig. 4). 
Hence, in tens of these m om en t arms and tire forces, the moment equations *re 


■ 1 X T* D 21 mW 009 T W + D 2 2 wW ain tW + D 2J^ 

- cos t w + D J2 p w sin t" + D |3 )]Ke« w (36) 


\ 
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w*** 


y 



CO*- 1 1 W 



k" « COORDINATES OF AXLE AT WHEEL w 
10 WHEN SHOCK STRUT IS IN FULLY 
EXTENDED POSITION 


R w = RADIUS OF WHEEL w 

Figure 4.- Orientation of landing gear relative to aircraft body axes. 

* t X T^31 ljW COS tW + ®32 wW 8 * n ^ + D 33 ) 

- xY(D xl u w cos t w + D 12 y w sin t w + D 13 )]K*6 W (37) 

R^ * l x Y(®nl jW cos T W + D 12 jj W sin t W + D l3 ) 

- Xj(D 21 y W cos T w + D 22 »i W sin t w + D 23 )]K*6 W (38) 

The forces relative to aircraft body axes are 



It should be noticed that, when the aircraft rebounds, 6 W - 0 and the 
friction forces and moments vanish. 

Static friction forces- Subsequent to wheel contact, wheel w is exposed 
to a static friction force. If this force is not sufficient to prevent slid- 
ing, a kinetic friction force will be generated. The kinetic friction force 
will persist until the velocity of wheel w in the plane of the landing plat- 
form is dissipated, at which point a static friction force is again encoun- 
tered. During static friction conditions, the velocity of the point of wheel 
contact *n the plane of the landing platform is zero. In this case, the wheel 
axle will move relative to the point of wheel contact by stretching the tire 
spring and the aircraft will oscillate under the restoring influence of the 
tire springs. 

If torsional stiffness is neglected, each tire may be treated as a three- 
spring system: a lateral spring that acts to restore deformations of the tire 

normal to the plane of the wheel, a spring that resists displacements in the 
direction of the strut axis, and a spring normal to these two directions. 

That is, if n is a unit vector normal to the plane of the wheel and s a 
unit vector coaxial with the shock strut, then the tire resists motions in the 
directions of n, s, and ft * s. 

In order to determine the spring forces, it is necessary to know the 
velocity components of the wheel axle in the directions of n, s, and n * s, 

n - l£ + MS + M£ 

s ■ ii + mj + nfe (39a) 

n * s = (nM - mN)i + (£N - nL)j + (mL - £M)k 

and the direction cosines £ w ,m w ,n w ,L w ,M w ,N w must be ascertained from the 
geometry of a particular aircraft. The tire displacements are 6 n normal to 
the plane of the wheel, 6y in the direction of the strut axis, and 6p normal 
to these two directions. Before touch-down, the displacement of the wneel axle 
relative to the point of wheel contact is zero. After touch-down, the velocity 
components of the wheel axle are 
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and 


w 

u a ’ 

/ ,W*Wv 

(u - l 8 ) i 

w 

v m 

a 

/ W.Wx 

(v - m s ) ; 

w 

w 

a 

w 

X ia 

<*10 - ^ * 

w 

X 2a “ 

, W w W N 

(x 20 -ms); 

< 

3a 


W.Wv 


, W W W* 

- (x 30 -ns) 


The x w io are the coordinates of the axle when the shock strut is In the fully 
extended position and s w is the shock-strut deflection at wheel w. 


I 

*. 


The forces generated by the tire displacements are 



n 

0 



where 


, w . w, w ,»> 

k “ k <p .<$ ) 

n n 


, v . w, w f w N 
k * k (p ,6 ) ; 

y y r ' 


! w . w, w *w. 
k - k (p ,6 ) 
P P 


w w 

and p and 6 are the pressure and tire deflection, respectively, at 
wheel w. 


,w 


If after touch-down wheel w has a velocity component Vp in the direc- 
tion of unit vector p, then a kinetic friction force will be generated. The 
magnitude of the force will be 


-w w_w _ w 
F ■ y Fj cos t ; 


w w, w. 
U » P (t ) 


r / i 


w w 

Likewise, a velocity component V n perpendicular to Vp will be accom- 
panied by a kinetic friction force F^ , where 

_w w_,w . w w w, w. 

F - y F, sin t ; y - y (t ) 
n 9 

Generally, these two components will exist simultaneously and give rise 
to the resultant kinetic friction force: 

Ff - y w F* 

The tire force components (F*j) relative to the landing platform are obtained 
by transforming the tire forces from strut axes to platform axes. These are: 
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A static friction force will persist until the friction force exceeds the 
spring force, that is, until 


(Fp 2 > (<F? t ) 2 + (F^) 2 ] 

When this condition is satisfied, the tire starts to slide and kinetic friction 
forces begin to operate. Relative to aircraft axis, the tire force components 
are 



When the aircraft touches down, that is, when the wheel contact condition 
(eq. (24)) is satisfied, the tire deflection rates are calculated and inte- 
grated and the corresponding spring forces determined. These are then com- 
pared with the forces required to produce sliding. When the kinetic friction 
forces exceed the spring forces, sliding ensues and the forces are modified 
accordingly. 

Effective tire stiffness- The effective tire stiffness at any orientation 
is a function of the component stiffnesses of the tire, the direction cosines 
defining the orientation of the shock strut and wheel relative to aircraft 
body axes, and the spatial orientation of the aircraft relative to the landing 
platform. The tire force produced by a tire deflection £ normal to the land- 
ing platform is K e 6R, where is the effective tire stiffness and ft is a 

unit vector normal to the landing platform. The corresponding tire ureas in 
the directions of n, s, and h * a are, respectively: 

(LD lt + MD 2J + KD || )«V[(p W ,6 W ) 

(*Dj j ♦ mDj , + nD #| )« w lJ(p w .« w ) 

KnM - mlODj, + (Ut - nL)D t> + (mL - MI)D || ]« ,, kJ(pV) 
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Therefore 
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•a 


-J 
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ak 
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K e ■ {<^13 + 3 + NDjsJX + ^13 + ® D 2 3 + n2 3J ) 2 Kj 

+ [(nM - mN)D 13 + UN - nL)D 2J + (mL - fM)D 33 ] 2 K^} l/z 

KINEMATICS AND DYNAMICS OF LANDING GEAR 


Landing Gear Forces 

Resultant force- The friction force and the platform reaction force are 
reacted at the wheel axle by a resultant force, which can be resolved into a 
component coaxial with the shock strut and a component normal to the shock 
strut. 

-y 

Axle force- The force F a acting at the wheel axle and the force Fg 
acting at the point of wheel contact with the landing platform combine with 
the gravity force M w g to Impart an acceleration to the wheel, that is, 

+ M v g - M"? 1 

where g is the gravity acceleration vector and r is the acceleration of 
wheel w (Mllwitzky and Cook, 1953). Therefore, 

p« - m w (? W - g) - (40) 

V 

The resultant force acting at the axle of wheel w has components f^ a 
relative to aircraft body axes. The force components normal to the shock 
strut may be expressed in terms of these components as follows. 

The scalar magnitude of the body axes components f ^ normal to the 
shock strut are 

f* a sin (cos -1 t w ) ; f* a sin (cos" 1 m w ) ; f* a sin(cos" 1 n w ) (41) 

In terms of the aircraft body axes components f^, the normal force FjJ 

has vector components as follows (Wills, 1958) a j 

fJJ - f" a (s * 1) * » + fj a (s X J) * s + f" a (s X fc) X i (42) j 

* * j 

Using vector identities, this equation can be rewritten as follows: 

fJJ - it\[i - (i* S)S] + f" a [3 - (j • S)S] + f, a lk - (k« S)S]} (43) i 

ft 

By substitution for § from equation (39a), equation (43) assumes the follow- 
ing form: 
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( 44 ) 


~ “ Amf 2a ” + l’ m * ,f Ta + ( l “ «*n)f* a ” ®nf* a ]j 

+ [-nif" a - nmf^ + (1 - nn)f* a ]k} 


Shock strut and platform forces- In addition to the force Fjj a normal to 

the shock strut and the axial strut force applied at the axle, the wheel 
is subjected to the platform reaction and friction force F^. This force is 

applied at the point of wheel contact and has components f* g relative to 
aircraft body axes. The resultant of these forces is 


F W - (Fj a + F^ + Fj) 

- {[(1 - M)f” a - lmfj a - *nf" a + lF a + f" g ]i + [-mlf? a + (1 - mm)fj a 
- mnf? a + oF^ + f? g ]j + [-ntf” a - nmf^ fl + (1 - nn)f” a + nF* + fj]k} 


V .**w 

- M (T 


g) 


This equation may be written more concisely as 


[F W ] - (L W ][F«] + F«| « W ) + IFJ] - M* ([X W ] - [g]> 



(45) 


(46) 
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From equation (40) , 


[F^J - -[Fp + M W ([X W ] - [g]) (47) 

Substituting from equation (39) yields 

[f£] - k"6 W [D][f”] + M W ([X W ] - [g]) (48) 



Substituting from equations (48) and (39) in equation (46) gives the equation 
of motion of the wheel mass tP: 



(49) 


Equation (49) gives the forces acting at wheel w and establishes a relation- 
ship between the forces due to tire deflection and the shock-strut forces. 
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The corresponding force acting on the aircraft is [pj] p where 


< F I>- 


+M*aiP'; - [,])> - 


,w /f ^w. 



( 50 ) 


and 



that is, the axial strut force Fg due to hydraulic resistance, air compres- 
sion, and internal bearing friction combines with and the prevailing 

aerodynamic, thrust , Inertia, and gravity forces to modify the aircraft motion 
subsequent to wheel contact. 


Influence of wheel mass- The fact that the mass of the wheel is a rela- 
tively small fraction of the total mass of the aircraft suggests that the 
wheel mass can be neglected without impairing the validity of the calculated 
results (Milwitzky and Cook, 1953). With this modification, the equations 
relating the tire forces to the strut forces assume the simpler form: 

/t w \ 

K*6 W lL w )lDHF*i ♦ F*| m W J - (5 

\ w / 

\° / 

By combining like terms, this equation can be rearranged as 


where 
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Matrix multiplication yields 


where 




*22 * “>12 
" ^11 


+ a D 


21 


+ nD 


22 


+ nD 
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+ nD 
+ nD 
4 nD 


ti 

32 

39 


Therefore. 


iy(« n u w CO. t w + ^ 2 u W .in T w ♦ *„) - - Fj(S.S) 


(52) 


For xero wheel ms.., the lending gear force, acting on the aircraft are 

/ 


pyiL W ]lD)[FjJ ♦ F; - I m 


Equation (51) permit, thi. equation to be written in the simpler form: 



(53) 


• llj) (54) 

Shock-strut closure and tire deflection rates- Subsequent to the time when 
*r " “< d n*xg + d t»*Tg + <*.»*?,> 

platform reaction forces, friction forces, and shock-strut forces ere gener- 
ated in accordance with the equation, formulated. Squat ion (32) gives the 


20 



relationship betveen the tire deflection forces and the '■hock-strut forces. 
The rate of tire deflection is & w . where 


;w . / w , w w . . . , w . w w . 

6 " d i>< u a + ^ x ja ' ^a* + d 2»< v a + "la * P x ja> 

+ d,,(w* + pXj a - qx^ a ) - (W + PXY • QX*) 


( 55 ) 


where 


and 


/ fl W*V N 

(u - l * ) 


vj - (v - m w . w ) ; 


W , W*Wv 

v a - <v - ti a ) 


x W (x^ 0 - i 


w . w ) ; 


2a 


(x 
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- .V) ; 


k ia 


/ w w w. 
(x J# -ns) 


Equations (52) and (55) can be solved to obtain ? w and s w , which can then 
be used to compute the platform reactions and friction fore i and the forces 
generated in the shock strut. When these forces are added to the aerodynamic , 
inertia, thrust, and gravity forces, the response of the aircraft subsequent to 
wheel contact can be determined. 


CONCLUSIONS 


The equations describing the friction forces generated during landing 
operations on ships at sea have been formulated. To simplify the formulation, 
it has been assumed that the force-deflection characteristics of the tire are 
linear, and that the behavior of the shock strut is relatively insensitive to 
variations in the air compression process. By ignoring variations in the 
polytropic exponent and by including only those characteristics of the tire 
and shock strut that contribute significantly to the generation of landing 
tear forces, a set of relatively simple equations is obtained. Nevertheless , 
these equations are considered adequate for practical purposes. 


t 


* 
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APPENDIX A 


TRANSFORMATIONS 


Transformation of Motion Vector Components 

A set of vector components in a coordinate system that is rotationally 
fixed is related to the components in the aircraft body axes by a transforma- 
tion equation of the form 


[A] - [TI^IEJ 


(Al) 


where 

[A] column vector of motion components in aircraft reference system 

[TJgA matrix that effects a transformation from fixed axes to aircraft body 
axes 

[E] column vector of motion components in fixed re c erence system 

Likewise, the components of a vector in the fixed reference system are 
related to the components in the moving ship reference system by a transforma- 
tion of the same form. That is. 


[SI - [T] es [E] 


(A2) 


where 

[S] column vector of motion vector components relative to ship axes 

[T] g£ matrix that effects a transformation from fixed axes to moving ship 

axes 

Similarly, a triad of ship axes components can be transformed to aircraft 
body axes by the transformation equation 

[A] - [T] sa [S] (A3) 


where 

[Tjg. matrix that effects a transformation from ship axes to aircraft body 

axes 

Substituting from equation (A2) in equation (A3) gives a transformation 
from fixed axes to ship axes, followed by a transformation from ship axes to 
aircraft axee: 


72 


(A4> 


1*1 ' I t 1sa it W e ' 

Finally, substituting from equation (Al) in equation (A4) yields the following 
matrix equation: 

[A] - [T W T W Ti ii A 

Therefore, ■ [I], where [I] is the unit matrix. Solving this 

matrix equation for [Tlg^ yields 

[T]g A ■ (A5) 

Since only orthogonal transformations are being considered, the inverse of a 
transformation matrix equals the transpose of the matrix and equation (A5) 
simplifies accordingly. That is. 


[T] 


-1 

ES 



(A6) 


where superscript T denotes transposition. 


Substituting from equation (A6) in equation (A5) yields the required 
transformation from ship axes to aircraft body axes. That is. 


[T] 


SA 


t T W T ^s 


(A7) 


In terms of the Euler angles v, 6, and 4 and with the conventional 
aeronautical rotation sequence, the required transformation matrices are 
(McRuer et al. , 1973) 



cos 6 cos 4 


cos 6 sin 4 


-sin 6 


IT] 


EA 


sin 4 sin 6 cos 4 
- sin 4 cos 4 


sin 4 sin 6 sin 4 
+ cos 4 cos 4 


sin 4 cos 8 


Icos 4 cos 4 sin 8 
L + sin 4 sin 4 


sin 4 cos 4 sin 8 
- cos 4 sin 4 


cos 4 cos 8 


(A8) 


For the transformation from fixed axes to ship axes, the Euler angles 
will be denoted by the capital Greek letters T, <3, and 4. In terms of this 
notation, the transformation matrix [T]gg asstaes a form identical to equa- 
tion (A8) 
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cos 0 cos ¥ 


cos 0 sin ¥ 


-sin 0 


r 


< 

t 


H 


4k 


IT] 


ES 


sin £ sin 0 cos f 
- sin ¥ cos $ 

cos ¥ cos $ sin 0 
+ sin sin $ 


sin V sin 0 sin £ 
+ cos V cos $> 


sin $ cos 0 


sin Y cos £ sin 0 cos 4> cos 0 
- cos Y sin 0 


(A9) 


The following equation gives the transposed form of this matrix: 


[T] 


ES 


cos 0 cos ¥ sin $ sin 0 cos V C 03 ¥ cos t sin 0 


cos 0 sin V 


- sin 0 


- sin V cos 9 

sin V sin 0 sin $ 
+ cos ¥ cos ♦ 

sin $ cos 0 


+ sin ¥ sin 4> 

sin ¥ cos $ sin 0 
- cos ¥ sin 4> 

cos ♦ cos 0 J 


(A10) 


A computer program, which solves the equations of the mathematical model 
of the aircraft, evaluates angles <|>, 6, and $ and hence determines the atti- 
tude of the aircraft as a function of time. The Euler angles are then used to 
compute the elements of the transformation matrix IT IgA* 

To determine the attitude of the ship, the components of the ship's 
angular velocity vector are measured and used to formulate the equations 
(McRuer et al., 1973) 

v 

P = $ - ¥ sin 0 


Q ■ 0 cos ♦ + ¥ cos 0 sin ♦ 

R * ¥ cos 0 cos 4> - 0 sin 4» 

• « • 

Solving these equations for 4, 0, and ¥, we obtain 
♦ » P + (Q sin & + R cos $)tan 0 
0 » Q cos ♦ - R sin ♦ 

¥ • (Q sin $ + R cos $)sec 0 


(All) 


(A12) 


The solution of these equations yields the required Euler angles <J> , 6, 
and ¥, which are then used to determine the elements of the transformation 
matrix [T]gg. After the formulation and transposition of this matrix, the 
product of [T]gA and [T]gg is formed. 

The transformation of vector components from ship axes to aircraft body 
axes is given by equation (A7): 

m SA • < T W T 4 


F / 


* 1 
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In terns of the direction cosines djj, relating aircraft body axes to 

ship axes* this matrix equation assumes the form 

'’Is*- 

The direction cosines d^j have the following values: 

d n “ ( cos 0 cos ^ c ° 8 0 cos ' + 008 0 8 * n ^ cos 0 8 *® * + 8in 0 sin 

d i 2 " f c ° 8 0 coa * sin ® co8 * “ sin ¥ cos *) 

+ cos 6 sin 4 (sin ¥ sin 6 sin 4 + cos ¥ cos 4) - sin 0 sin 4 cos 0] 

d J3 » [cos 6 cos 4(cos ¥ cos 4 sin 6 + sin ¥ sin 0) 

+ cos 8 sin 4 (sin ¥ cos 4 sin 0 - cos Y sin 4) - sin 0 cos 4 cos . 0 ] 

d 2 - [ (sin $ sin 0 cos 4 - sin 4 cos 4)cos 0 cos ¥ 

+ (sin 4 sin 0 sin 4> + cos 4 cos 4)cos 0 sin Y - sin 4 cos 0 sin 0 ] 

d 22 ■ [(sin 4 sin 0 cos 4 - sin 4 cos 4) (sin * sin 0 cos ¥ - sin ¥ cos 4) 

+ (sin 4* sin 0 sin 4 + cos 4 cos 4) (sin ¥ sin 0 sin 4 + cos ¥ cos 4) 
+ sin 4 cos 0 sin 4 cos 0] 

d 2S ■ [(sin 4 sin 0 cos 4 - sin 4 cos 4) (coa ¥ cos 4 sin 0 + sin ¥ sin 4) 

+ (sin 4 sin 0 sin 4 -I- cos 4 cos 4)(*in Y coa 4 sin 0 - cos ¥ sin 4) 
+ sin 4 cos 0 cos 4 coa 0] 

d sl “ [(cos 4 cos 4 sin 0 + sin 4 sin 4) (cos 0 cos ¥) 

+ (sin 4 cos 4 sin 6 - cos 4 sin 4)cos 0 sin ¥ - cos 4 cos 0 sin 0 ] 

d j 2 ■ [(cos 4 cos 4 sin 0 ♦ sin 4 sin 4)(sln 4 sin 0 cos ¥ - sin ¥ cos 4) 

♦ (sin 4 cos 4 sin 0 - cos 4 sin t)(aln ¥ sin 0 sin 4 + cos ¥ cos' 4) 
+ cos 4 cos 6 sin 4 cos 0] 


d d d 

111213 


d d d 
2 1 22 2 8 


d 31 d S2 d 3 3 


- ID] 
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d 3 3 ■ l (cos \p cos <fr sin 0 -*■ sin ip sin ^(cos Y cos * sin 0 + sin Y sin #) 

+ (sin i|> cos <p sin 0 - cos ip sin $)(sin Y cos * sin 0 - cos Y sin *) 

+ cos Ip cos 0 cos $ cos 0] 

Although the direction cosines are useful for transformation purposes, 
they are not convenient measures of aircraft attitude. A conversion from 
direction cosines to a set of Euler angles that represents the attitude of the 
aircraft relative to the ship can be effected by the method described by 
Meyer et al. (1967). 


For the conventional aeronautical rotation sequence, a rotation matri;: 
can be generated ae the product of three rotation matrices as follows: 


where 


(T 


[D] - tT 1 (0 1 )][T 2 (0 2 )]![T,(0 3 )] 


r ° 

1 (0 1 )] ■ I 0 cos 0 X 
\0 -sin 0 X 



cos 0 


(A13) 


11 ,( 6 ,)) 


cos 0, 0 -sin 0 


1 

0 




( cos 
-sin 


cos 0 3 sin 0 3 

0 , cos 0 . 


0 

, 0 

0 1 

The product matrix, equation (A13), yields the following direction cosine: 
d lx - cos 0 2 cos 0 3 
d 12 ■ cos 0 2 sin 0, 
d xs - -sin 6 2 

d 2X » sin fl x sin 0 2 cos 0, - sin 0 3 cos 0 X 
d 22 ■ sin 0 3 sin 0 2 sin 0 X + cos 0 3 cos 0 X 
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sin 0 X cos 0 2 





> 
v > 



I 


/I 



f 


d 31 ■ cos e 3 cos e x sin e 2 + sin 0 3 sin 0 X 
d 32 - sin 0 3 cos 0 X sin 0 2 - cos 0 3 sin 0 X 
d 33 - cos 0 X cos 0 2 


The following combinations of these equations are required to convert the 
direction cosines d^j to the Euler angles 0 X , 0 2 , and 0 3 : 


d xx sin 0 3 - d X2 cos 0 3 » 0 (A14) 
d 31 sin 0 3 - d J2 cos 0 3 ■ sin 0 X (A15) 
d 2X sin 0 3 - d 22 cos 0 3 ■ -cos 0 X (A16) 

d X3 - -sin 0 2 (A17) 
d 23 sin 0 X + d J3 cos 0 X ■ cos 0 2 (A18) 


From equation (A14), 


tan 0 3 



(A19) 


Equations (A15) and (A16) give 


tan 0 


or 


tan 0 


[ d 3x sin 0 3 - d 32 cos 0 3 "j 
d 22 cos 0 2 - d 2x sin 0 3 J 

f ~ d 3 1 tSn 9 3 - d 3 2 1 

1 [_ d 22 d 21 tan ®*J 


From equations (A17) and (A18) 


„ d ll 

tan 0, ■ - -x 3 — t — r—r — 

2 d 23 sin 0 X + d 33 cos 6 X 


Given the nine direction cosines and the rotational sequence, the three Euler 
angles 6 X , 0 2 , and 0, can be computed. For the present application, 0, 
corresponds to the aircraft yaw angle relative to the landing platform. 

Anglss 0 2 and 0 t are pitch and roll angles, respectively, relative to the 
landing platform. It should be noted that the computed values of 0^ are not 
unique since tan 0 is a many valued function, that is. 
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tan 0 ■ tan(nir + 9 ) 


where n la a positive or negative integer. However, for the case being 
considered and for the angles anticipated, only these solutions corresponding 
to n » 0 will be required. 


APPENDIX B 


INTERNAL SHOCK-STRUT FORCES 


Hydraulic Forces 

The hydraulic force F^ is obtained by making use of the equation for 
the discharge through an orifice (Mllwitzky and Cook, 1953) 

0 • C d*„Vf <»h ' Pa> < Bl > 

Q volumetric rate of discharge 
C d coefficient of discharge 
A n net orifice area 

hydraulic pressure in lower chamber 
p air pressure in upper chamber 
p mass density of hydraulic fluid 

Using continuity considerations, the volumetric rate of discharge can be 
expressed in the alternative form: 


Q-AhS 


(B2) 


where A}, is the hydraulic area and S is the telescoping velocity of the 
shock strut. Equating these two expressions for the volumetric rate of dis- 
charge yields the pressure differential 


P h 


p . 


PA h S 

2 < c dV a 


(B3) 


The hydraulic resistance due to the telescoping of the strut is given by 
the product 


pA? 

<Pw - P.)Av S 2 - F. (B4 

^ 2(c dV h 

This equation can be used for both the compression and elongation strokes by 
Introducing the factor S/|£| to Indicate the sign of the hydraulic force as 
follows: 
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(B5) 



PA h 

2(C dV 


S 2 


Pneumatic Forces 

The pneumatic force F a in the upper chamber depends on the inflation 
pressure and the exposed area. It is assumed that the air pressure obeys the 
generalized gas law (Milwitzky and Cook, 1953) 

pv n » constant ■ p Q v^ 


or 


P 



where 

p air pressure in upper chamber of shock strut 

p Q air pressure in upper chamber for fully extended strut 

v air volume of shock strut 

v 0 air volume for fully extended strut 


(B6) 


The instantaneous volume v is equal to the initial volume minus the swept 
volume, where the swept volume is the product of the pneumatic area A a and 
the stroke of the strut S, that is. 


v - (v 0 - A a S) 

Substituting this value in equation (B6) gives 

- • 4rM’ 

It follows that the pneumatic force Is 

P a " P^a " PjA 




(B7) 


(B8) 


The exponent n depends on the rate of compression and rate of heat transfer 
from the air to the surrounding environment. Low rates of compression corre- 
spond approximately to the Isothermal case and a value of n ■ 1 . Higher 
rates of compression approach the adiabatic condition and a limiting value of 
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n » 1.4. A limited amount of experimental data obtained in drop teats indi- 
cates that, in most practical cases, a value of n » 1.1 may be used 
(Milwitzky and Cook, 1953). 


Internal Friction Forces 

According to the law of friction, the coefficient oi friction u is 
defined as the ratio of the friction force to the normal force. This coeffi- 
cient is somewhat greater under conditions of rest (static friction) than under 
conditions of sliding (kinetic friction). Machine designers usually classify 
frictional resistance as friction between dry surface, friction between imper- 
fectly lubricated surfaces, and friction between perfectly lubricated surfaces. 
The internal friction in landing gear shock struts usually Involves relatively 
high normal pressures and small sliding velocities. Moreover, the usual types 
of hydraulic fluid used in shock struts have Imperfect lubricating properties. 
Therefore, for the sake of completeness, it will be assumed that the Internal 
friction between the bearings and cylinder walls approaches the dry friction 
condition. Hence, the internal friction forces which depend on the magnitude 
of the bearing forces, the orientation of the gear, the spacing of the bear- 
ings, and the appropriate coefficients of friction are obtained as follows 
(fig. 5). 




The bearing forces are related to the force FjJ a normal to the axis of 

the shock strut. This force Is required to determine the Internal friction 
forces ff, where 


f 


w 

f 



(B9) 


where 

f f axial friction force for strut at wheel w 

u? coefficient of friction for upper bearing (attached to inner cylinder) 

F* normal force on upper bearing (attached to inner cylinder) 

coefficient of friction for lower bearing (attached to outer cylinder) 
normal force on lower bearing (attached to outer cylinder) 
sign of friction force 




_S_ 

.Si 


(ij - s) 

F W - f W — 

1 N a (t« ♦ S) 

A “ + F T> - 


s |p‘“ - s > , . * .1 

■ w F «j(7 *77> <u ‘ " ,> u j 


(BIO) 


(see fig. 5). 


Total Strut Force 

The sum of the pneumatic* hydraulic* end lnterral friction forces gives 
the total shock-strut 'force* which la 

’• • [ft 5^7 *' * 4V * ifi lllfri) <* • •*> * - ; |] 


(Bll) 
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